Cosmic microwave background bispectrum on small angular scales 
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This article investigates the non-linear evolution of cosmological perturbations on sub-Hubble 
scales in order to evaluate the unavoidable deviations from Gaussianity that arise from the non- 
linear dynamics. It shows that the dominant contribution to modes coupling in the cosmic microwave 
background temperature anisotropics on small angular scales is driven by the sub-Hubble non- 
linear evolution of the dark matter component. The perturbation equations, involving in particular 
the first moments of the Boltzmann equation for photons, are integrated up to second order in 
perturbations. An analytical analysis of the solutions gives a physical understanding of the result 
as well as an estimation of its order of magnitude. This allows to quantify the expected deviation 
from Gaussianity of the cosmic microwave background temperature anisotropy and, in particular, to 
compute its bispectrum on small angular scales. Restricting to equilateral configurations, we show 
that the non-linear evolution accounts for a contribution that would be equivalent to a constant 
primordial non- Gaussianity of order /^^ ~ 25 on scales ranging approximately from i ~ 1000 to 
£ ~ 3000. 

PACS numbers: 98.80.-k 
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I. INTRODUCTION 



The cosmic microwave background (CMB) offers a 
unique window on the physics of the early Universe, and 
in particular on inflationary models. The angular power 
spectrum of the CMB anisotropics has been extensively 
used to set constraints on the shape of the inflationary 
potentials; see e.g. Ref. [i|. The statistical properties 
of the temperature anisotropies and polarisation depend 
both on the inflationary period during which they were 
created and on the physics at play after Hubble-radius 
crossing and during the recombination. At linear order 
in metric perturbations, those latter physical processes 
amount to affect the metric perturbations by a multi- 
plicative transfer function. The characteristic features 
observed in the temperature anisotropy spectrum origi- 
nate from the development of acoustic oscillations that 
this transfer function encodes. The overall amplitude 
of the metric perturbation and its scale dependence are 
however determined by the inflationary phase. 

At linear order, the calculation of the transfer function 
- and hence the detailed shape of the temperature power 
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spectra - for generic inflationary models requires the iden- 
tification of the relevant degrees of freedom during infia- 
tion (see e.g. Refs. [2, 3, 4]), as well as a full resolution 
of the dynamics up to recombination time. All these as- 
pects are now fully understood (see e.g. Refs. 0, '6'| and 
references therein). 

At this level of description, the metric perturbations 
are linearised so that the non-linear couplings that are 
inherently present in the Einstein equations are ignored. 
Therefore models that predict Gaussian initial metric 
fluctuations are expected to induce cosmic flelds with 
Gaussian statistical properties. This is a priori the case 
for generic models of inflation. It has to be contrasted to 
models with active topological defects - such as cosmic 
strings - that have soon been re cog nised as a source of 
large non-Gaussianities 0, [a, [3, [Ifl]. The current data 
however clearly favor only mild non-Gaussianities al- 
though those might be larger than those induced by pure 
gravity couplings. This is not the case for single fleld 
slow-roll inflation for which it has been unambiguously 
shown in Ref. [ll| that it can produce only very weak 
non-Gaussian signals, that are bound to be overridden 
by the gravity induced couplings. It has however been 
realised that some models of inflation might produce 
significant deviation from Gaussianity in the context of 
multiple-field infiation p,|l3„ 14, 15,Jll [O, D,8, 19, 2^ 
or with non-standard kinetic terms [2l|. The question 
of the observation of primordial non-Gaussianities is 



largely open. 

In general however, primordial deviations from Gaus- 
sianity are in competition with the couplings induced 
during the non-linear evolution of the cosmic fields. It 
has triggered general studies aiming at characterising the 
bispectrum to be expected in the observation of the cos- 
mic microwave background temperature anisotropics and 
polarisations whether it arises from inflation or from sub- 
sequent effects. 

This task is multi-fold. It requires a proper identifica- 
tion of the mode couplings (at the quantum level) during 
the inflationary phase - so taking into account the usual 
gauge freedom - as well as a second order treatment of the 



post-inflationary evolution. While the former has been 
set on firm ground [ll|, [22] , the latter issue is still largely 
unexplored. This article proposes both numerical and 
analytical insights into it. 

Hereafter, we assume that on super-Hubble scales, the 
only significant scalar perturbations are adiabatic and 
that they obey a nearly Gaussian statistics [43|. To be 
more precise, they are described in Fourier space by a 
single variable Co(k), k being a comoving wave- number, 
that satisfies 



(Co(ki)Co(k2))=(5D(ki+k2)Pc(^i) 



(1) 



and 



J 



(Co(ki)Co(k2)Co(k3)) = 2fo(ki +k2 +k3) /«Jki,k2)Pc(fci)n(fc2) +sym. 



(2) 



where "sym." stands for the two other terms obtained by permutation of the wave-numbers. This defines the primordial 
power spectrum P|^(fc) and the primordial mode coupling amplitude [48| f' . Considering an observable quantity 9 
related to the perturbation variables, the effect of evolution can generically be recapped 49] as 



e{k) = Tj'\k) Co(k) + J "^'^^P' fe(k-ki-k2)Tf '(ki, k2) Co(kl)Co(k2) 



(3) 



where Tg (fc) is the linear transfer function and Tg (ki, k2) is the second order transfer function. 9 can be thought 
as being e.g. the observed temperature anisotropics, but it could also stand for the CMB polarisation, or even cosmic 
shear surveys. When computing the bispectrum of 9 there will be a contribution from the mode couplings induced 
by the second order transfer function T'^' and the possible initial non-Gaussianities, 



(0(ki)0(k2)0(k3)) = 2fc(ki -f k2 +k3) [4^^(ki,k2) +/;^,(ki,k2)] Tg^'^{h)Tg''Hk2)V'\ks) Pdh)Pc{k2) 



-(1)/ 



-(1)/ 



-|-sym. 



(4) 



r 



where f ^ is related to the second order transfer function 
by 



-(2) 



r;^^(ki,k2)^/: fki,k2)r;^^(|ki + k2|) 



(1)/ 



(5) 



-(2) 



The full derivation of the details of Tg is a fantas- 
tic task. It requires an understanding of the metric 
fluctuations behaviour at second order, from radiation 
dominated super-Hubble scales to matter dominated era 
at sub-Hubble scale, as well as a comprehension of the 
physics of recombination - through the Boltzmann equa- 
tion - at a similar order. Such a task has been under- 
taken by several authors 50] and the multitude of effects 
at play needs to be sorted out. So, the goal of this article 

(2) 

is not to provide an end to end calculation of 7^ , but 
to show that on small scales one can extract the domi- 
nant terms in order to get an insight into this physics at 
second order. 

Modes coupling due to gravitational clustering is, by 
far, not a novel subject. It can be traced back to early 
works by Peebles [23| where the function T^^^ for the 



non-linear sub-Hubble evolution of cold dark matter field 
(CDM) during a matter dominated era was derived. Gen- 
eral modes coupling effects, within the same regime, 
has been extensively studied in the eighties and nineties 
where a whole corpus of results has been obtained (see 
e.g. Ref. [2^ for an exhaustive review). On sub-Hubble 
scales, the second order mode coupling function for the 
gravitational potential reads 



/*Jki,k2) = 






_iJ2a2|ki- 
1 ki • k2 , 1 ki 
2 



-k2| 
k2 



fc2 



^2 



2(ki 

7 



''■1"'2 



(6) 



in the particular case of an Einstcin-dc Sitter universe 
(here a is the scale factor and H the Hubble parame- 
ter). This well-established result proved to be useful for 
observational cosmology. The angular modulation it ex- 
hibits has indeed been observed in actual data sets; see 
e.g. Ref. |25|. 

The fact that on sub-Hubble scales, that is k"^ ^ 



H^a^, the non-Gaussianity is driven by the non- 
hnearities of the CDM sector, and that they can start 
developing even before equahty is one of the leading ideas 
of the present study. Indeed temperature anisotropies on 
small angular scales ~ i.e. beyond the first acoustic peak 
- mostly trace the gravitational potential [5i| long after 
it has entered a sub-Hubble evolution and already during 
the matter dominated era. It is then natural to expect 
that the temperature anisotropies should be substantially 
determined by a form close to that of Eq. ^ . 

The goal of this paper is to evaluate how close we are 
from the behaviour ^ depending on scales, to which 
extent the temperature anisotropies trace this form and 
finally to estimate the amplitude of the temperature bis- 
pectrum on small angular scales. In this work two ap- 
proaches will be compared: a full numerical integration 
of the second-order equations presented in § |lll where 
the main approximation lies in the modelisation of the 
Compton scattering collision term at second order, see 
Eq. (|5D|) . and an approximate analytical resolution dis- 
cussed in § mil 

The bottom line of our analysis is that on small an- 
gular scales, the density perturbation of the cold dark 
matter starts to dominate the Poisson equation so at the 
time of decoupling we can assume that the system is split 
in (1) the evolution of CDM and (2) the evolution of the 
photons-baryons plasma which develops acoustic oscil- 
lations in the gravitational potential determined by the 
CDM component. As we shall demonstrate, at second 
order the dominant term of the temperature fluctuations 
is driven by the second order gravitational potential. Our 
approximation requires to consider a regime in which the 
Silk damping is efficient, that is wave-modes larger than 
the damping scales, hence corresponding to multipoles 
roughly larger that 2000. This picture will be shown 
to be in agreement with the numerical estimation (see 
§ nil Dp . We then proceed in § IIVI by a computation of 
the bispectrum in which wc show that, for equilateral 
configurations, the non- linear dynamics has an ampli- 
tude equivalent to that of a primordial non-Gaussianity 
with constant /^j^ of order 25. A back-of-the-envelop ar- 
gument allows us to understand the magnitude of this 
number. 



A. The background dynamics 

The background space-time is described by a 
Friedmann-Lemaitre metric with scale factor a and cos- 
mic time t. It is convenient to rescale the scale factor 
such that 

y = pnjpr , 

where pm and pr sue the background matter and radi- 
ation energy densities respectively. The matter energy 
density can be decomposed as the sum of a cold dark 
matter component, that does not interact with normal 
matter, and a baryonic component, that can be coupled 
to radiation by Compton scattering prior to decoupling. 
We thus set Pm = Pc + Ph where pc and ph refers to CDM 
and baryons respectively. It follows that 



Pin = 



1-/. 



with /j^ = ribo/^mo ~ 0.18. The Friedmann equation 
then takes the simple form 



n — /toq 



i + y 

22/2 



(7) 



when we neglect the contributions of the spatial curva- 
ture and of the cosmological constant, which are negli- 
gible for the whole history of the Universe but very re- 
cently. Ti = a! I a is the conformal Hubble parameter and 
a prime refers to a derivative with respect to the confor- 
mal time T] defined by dt = aAT\. Ticq is the value of the 
H. at equality, that is when y = 1. 

The equation of state of the background fluid, com- 
posed of a mixture of non-relativistic matter and radia- 
tion, is w — 1/[3(1 -t- y)] and the density parameters of 
matter and radiation are 



i ^m 



y 



1 + 2/ 



n, = 



1 



l + y 



(8) 



and indeed ilc — ^m(l ~ A)- 

Equality takes place at y = 1, from which we deduce 
that 



j/o = 1 + Zcq = 3612 6^7 ' 



0.15 / ' 



(9) 



II. PERTURBATION THEORY 



This section is devoted to the presentation of the per- 
turbation equations, up to second order, and of the initial 
conditions used in our study. We set the main notation 
and describe the background dynamics in § III Al and we 
define the perturbation variables in iJ IIIBI The perturba- 
tion equations and initial conditions are then presented 
in § III CI and IIIDI respectively. 



where 62.7 = To/2.7 K is the temperature of the CMB 
today, and h the value of the Hubble constant in units 
of 100 km/s/Mpc. Equation ([7]) evaluated today implies 
that Tieq ~ 7io-\/22/o so that 

T^cq -- 0.072 ano/i^ Mpc"^ (10) 

The last scattering surface corresponds to a redshift 1] 

l + ^Lss = 1090±l = j/o/yLSs, (11) 

and is mildly dependent of f2cO £ind rJbo- This implies 
that ?/LSS ^ 3.3. 



In Fourier space, a mode is super-Hubble when krj <t^ 
1 and sub-Hubble otherwise. The mode becoming sub- 
Hubble at equahty corresponds to a comoving wavelength 
of 






if we choose units such that oq = 1. 
We also introduce the parameter 



Mpc 



(12) 



R 



3pb 

4pr 



^/b y > 



(13) 



which will be useful to describe the physics of the 
baryons-photons plasma. 



B. Perturbation variables 

We focus on the dynamics of scalar perturbations (see 
e.g. Refs [2a, [23, [23| for analysis of vector and tensor 
modes generation at second order). In Newtonian gauge, 
we can expand the metric as 



ds' 



a^{T]) [-(1 + 2$)dry2 + (1 - 2*)<5yda;MxJl 



(14) 
where $ and ^ are the two Bardeen potentials. 

The various fluids contained in the universe will be de- 
scribed at the perturbation level by their density contrast 
6 and their velocity field. For the latter, we decompose 
the time-like tangent vector to the fluid worldlines ac- 
cording to 

a 

where the first term accounts for the background Hub- 
ble flow. The perturbation v'^ is further decomposed as 
v^^ = (ti°,ti*) with u* = d^v and w" is constrained by the 
normalisation Ufj_u'^ = — 1. 

When dealing with perturbations beyond first order, 
we assume that the perturbation variables are expanded 
according to 



1 



X = X(1) + -X 



(2) 



(15) 



where X'^^'> satisfies the first order field equations while 
the second order equations will involve purely second 
order terms, e.g. X^'^' as well as terms quadratic in 
the first-order variables, e.g. [X^^^]^. Thus, there shall 
never be any ambiguity about the order of perturbation 
variables involved as long as we know the order of the 
equation considered, and consequently we will usually 
omit the superscprit (1) or (2) which specifies the order 
of the perturbation. For a general discussions on sec- 
ond order perturbations and gauge issues, we forward to 
Refs. [2i,[33,[3l 



C. Perturbation equations 

In this article we shall focus on the CDM-radiation- 
baryons system. Each component has a constant equa- 
tion of state, that is P = wp with w' = so that c^ = w. 
Non-relativistic matter is described by a pressureless fluid 
with Wra = Wh = Wc = and radiation satisfies w^ = \. 

In full generality, the evolution of each component can 
be obtained from the Boltzmann equation satisfied by 
the distribution function fa{x^ ,pv) for this matter com- 
ponent. The stress-energy tensor can then be defined by 
integrating over momentum as 



where 7r+ [p) is the volume element on the tangent space 
in a;" such that p^ is non-spacelike and future directed 
(see e.g. Refs [Mill) 

The first moment of the Boltzmann equation then gives 
a conservation equation of the form 32] 



v^rr 



F" 



(16) 



where F^ describes the force acting on the fluid labelled 
by a and satisfy F^u^ = and ^^ F^ = 0, which is noth- 
ing but the action-reaction law (equivalently obtained 
from the Bianchi identity). Projecting along and perpen- 
dicular to u'^, we can extract respectively the continuity 
and Euler equations. 



1. Linear order 

Linear order calculations are used in particular to set 
the source terms of the second order equations. We 
closely follow the standard calculations and the main in- 
gredients are recalled here. At linear order, the continuity 
equation for a fluid labelled by a takes the form [S, [a [0| 

5'a + 3W(c' „ - wa)5a + (1 + Wa) (Az;, - 3^') = , (17) 



while the Euler equation 



nil-ici 



$ 



1 



^a - l^^a , (18) 

6 



where tt^ is the contribution of the anisotropic pressure. 
In deriving Eq. (fT^ . we have decomposed the force term 
as Fi = d^F = d.,[{p + P)T] and Fo = 0. Since F^ 
vanishes at the background level, Ta is gauge invariant. 
Wa and Cs.a are respectively the equation of state and 
sound speed of the component a. 

In our analysis, we consider three components. Dark 
matter (label c) is described by a perfect fluid with Wc = 
and TTc = interacting only through gravity (jFc — 
0). Baryons and photons are coupled through Compton 
scattering so that T^ and J> do not vanish. The action- 
reaction law (or equivalently the conservation of the total 



Component 


w 


J^ 


TV 


description 


CDM (c) 











fluid 


Baryons (b) 





n 





fluid 


Photons (r) 


1 

3 


J^r 


TTr 


kinetic 
(8 moments) 



The Einstein equations reduce to the set 



TABLE I: Summary of the properties and descriptions of the 
matter components considered in our analysis. 



stress-energy tensor of matter) imphes that F^. 
from which we deduce that 

J^i- ^ — R-Fh ■ 

At hnear order, it is easily shown that 

^r = t'{v\, - Vr) , 

where 

r = aueCTT , 



(20) 



(21) 



with Tie being the free electrons number density and ctt 
the Thomson scattering cross-section. It follows that 
baryons will be described by a fluid [wh = = TTb = 0) 
interacting with radiation. In general radiation enjoys a 
non-vanishing anisotropic pressure (tti- ^ 0) and should 
actually be describe by the full Boltzmann hierarchy. For 
the linear order calculations we choose to extend the 
fluid description by including the eight first moments of 
this hierarchy, including polarisation (see Appendix Elfor 
these equations). Our choices for the modclisation of the 
matter sector are summarised in table HI 



A* - 3H*' - iu^^ - -n^ J2 ^"'^^ = ^ ' (22) 



*" + n^<i> + -A($ - *) + n<^' + 2H*' 





+27i'$ - -n^nA- = , 




(23) 


-Fb, 


* - $ = n.n'^T:, , 




(24) 




^'+H^+'^H^ ^ na{l+Wa)Va-~ 


= . 


(25) 


(19) 


a—r.m 







When the anisotropic pressure can be neglected, and in 
particular in the tight coupling regime discussed below, 
Eq. ([24| implies that $ = \1/. Note that in this analysis 
we actually ignore the neutrinos the effect of which is 
thought to be marginal on the qualitative results we will 
obtain. 



2. Second order 

At second order, any first order equation, schemati- 
cally written as I?[Ar*^^^] — 0, of the first order perturba- 
tion variables X^^^ will take the general form 

x>[a:(2)] = s 

where S* is a source term quadratic in the first order vari- 
ables. 

For the continuity and Eulcr equations respectively 
read 



J 



Se,a - 2(1 + Wa) {6**' - ^AVa ~ d,Va [(1 " Sc^ JHS^^;, + 29'< + 29'$ - 39**] } 



1 + Wa 



Oibe a 



l+Wa 



[{Sad^VaY + n{l - iWa)5ad^Va + 5ad,^ + 27i(l - Sc^ J($ + 2^)d,Va 



'2<^^^v' + 4$(9,$ 



+ iWd.Va + 4*a,t;; - 29, {&Va^^Va) + 2^ 



(cL)' 






SadiVa 



(26) 



(27) 



As long as CDM is concerned, these source terms and Eqs. (|17m8l) gives the full second order evolution of the fluid. 
As already seen at first order, the fluid equations for the baryons and photons must include interaction terms, that 

(2) (2) 

is Tr and T^ , that derive from the Compton scattering collision term entering the Boltzmann equation for the 
radiation. 

In the baryon rest-frame, this collision term includes only two types of contributions [33, H, ISa] • First, there is a 
term involving first order perturbation quantities and whose form is 



2 \ Tie dT Xe 



(28) 



where C^^^ is the first order collision term. This contribution involves the fluctuation of the visibility function, 
that is of the electron density n,, and of the ionisation fraction x^. It accounts for the fact that a hotter or denser 



region decouples later. Its typical magnitude is of order iC^^'^Sne/ne 
perturbation variables and whose form is 



Second, there is a term involving second order 



2 L 



(29) 



The forces derived from these two terms will satisfy by construction the action-reaction law P^ . and this holds 
in any reference frame and at any order. This explains why the computation is easily carried out in the baryons 
rest-frame 34 1. 

Then, when changing frame from the baryons rest-frame to the cosmological frame, where the computations are 
actually carried out, a second series of terms appears. They are of the form r'/^*^^ x [w^^^]^ and r'f^^^ x [u'-^-'], where 
f^^' and f^ are the background and first order distribution functions as well as similar terms for the polarisation 
(see Ref. [34| for the exact form of these terms.). 

Now, the contribution (|28p is proportional to the collision term at first order. This implies that it will thus be 
negligible as long as tight coupling between baryons and photons is maintained at first order, i.e. as long as r'/fc 3> 1. 
We are thus left only with the contribution (j29p . This term thus enforces the tight coupling regime at second order. 
Then, as long as tight coupling is effective, it is obvious that the second series of terms arising from the change of 
frames should compensate each other to give a vanishing contribution. We shall thus model the interaction term 
entering the Euler equations by 



r 



(2)^_^(l)[;f(2)] 



(30) 



that is by assuming that it keeps the same functional form as at first order. In conclusion, the continuity and Euler 
equations for baryons and photons with the interaction term ([50]) are the exact fluid limit of the full Boltzmann 
equation at second order as long as tight coupling is effective. This implies that at second order, and similarly as at 
first order, the two tightly coupled fluids are equivalent to a single perfect fluid; see section UlI Al Again, this stems 
from the action-reaction law which implies that Eq. (J19p has to hold at any order and in any reference frame and that 
there cannot appear any external force acting on the resulting effective fluid since it is only coupled to other matter 
components through gravitation. 

When r'/fc becomes of order unity, tight coupling stops being effective. But thanks to Silk damping, the terms of 
the form r'/^°-' x [w^^^]^ and r'f'^^^ x [w*-^-'], will still be negligible compared to the one we kept to obtain Eq. (|5D|) . 

Let us now turn to the Einstein equations (|22ll25p . Their source terms read respectively 

Si = -8*A* - 3ai*a** - 3*'2 + 3^2 Y^ na{l+Wa)^^Va^'Va-l2n^■i'^ , (31) 

a— r.c.b 

5^2 = AH^'i'^ +'^di'i'd'^ + ^-^A-^ + 8n^-^' + 8n'-^^ + ^'^ +n^ ^ na{l + Wa)d^Vad^Va , (32) 



So = -m^ - A- 



25, *5'* + m^ 



-3{AA)-^d^dj 



fia(l + Wa)diVad''Va 



2d'^d^^ + 3n^ ^ na{l+Wa)d'Vad^Va 



S4 = 2n^'-i^-^' + 2d-\^'d,-^)+3n'' Y. ^adr^ [{1 + Wa)'i'd,Va ~ {1 + cljSad.V,] 



(33) 



(34) 



r 



This provides all the source terms appearing at second 
order. 



3. Note on our conventions 



the source terms can be written as a convolution on the 
wave-numbers ki and k2 such that ki -I- k2 = k. For 
simplicity of notations we will write the integral factor of 
the convolution as 



Since the first and second order equations are conve- 
niently solved in Fourier space, the quadratic terms in 



d^kid^ks ,. ^. ,, 

' (2^)3/2 fo(I^l+k2-k) 



We also define /i = ki •k2/(A:i/c2)- Unless explicitly spec- 
ified, we choose the convention of the Fourier transform 
in which the factors of (27r)"/^ are symmetric for the 
Fourier transform and its inverse, n being the dimension 
of space. 



D. Initial conditions 

To integrate this system of equations, we need to set 
the initial conditions both for the first and second or- 
der variables deep in the radiation era for super-Hubble 
modes at the initial time Tyinit, that is modes such that 
krjinit < 1- 



First order 



Deep in the radiation era, this implies that 

for modes such that kijinit ^ 1- Since the density con- 
trast of the total fluid 



Or + :r- <^m , 



1+y 1+y 



(37) 



6 ~ Si deep in the radiation era. From Eq. ([22]) we deduce 
that 

(5(i)(k,77i„it) = -2$(k,77i„it) . 

Now, assuming adiabatic initial perturbations, we must 
have 



At first order, we rely on the comoving curvature per- 
turbation which is constant on super-Hubble scales. It 
is well-known [a, 101 that for a perfect fluid with a time- 
dependent equation of state (c^ 7^ w), the comoving cur- 
vature perturbation, defined by 



7^(1) ^ $(1) 



3{i + w)n 



^'(1) ^n^ 



(1) 



(35) 



is conserved on super-Hubble scales for adiabatic pertur- 
bations. Inflationary models predict the initial power 
spectrum of T?,^^-* on those super-Hubble scales from 
which one can deduce the power spectrum of the gravita- 
tional potential. If we choose jyinit such that the decaying 
mode is negligible and \[''^^-' is constant on super-Hubble 
scales then, still neglecting the anisotropic pressure. 



7^(l)(k,r;init) = 



5 + 3w, 
3 + 3w 



*^'Hk,r7i„it) 



(36) 



5(i)(k,r?i„it) = 5i'^(k,ryinit) = -J^^^k, r/mit) ■ 
Since baryons and photons are tightly coupled deep in 
the radiation era, we deduce that 



= -2*(k,'7mit) • 

This completely fixes the initial conditions for the set of 
first order perturbation equations. 



2. Second order 

At second order the previous procedure can be gener- 
alised [3a, [33 . It was shown that on super-Hubble scales 
and for adiabatic perturbations, the variable 



n 



(2) ^ -^(2) 



3{l + w)H\ n 



[I + icl) 



3(1 + w) 



3(1 + w) 



(5* 



(38) 



is a conserved quantity on super-Hubble scales for adiabatic perturbations. Once the decaying modes are negligible 
so that <I>^^' and ^'^•' are constant, we can express them in terms of HS'^^ as 



*(2) ^ ^^i3(l + ^„)7^(2)+44±^vI;(l)2_2A-l 
5 + 3w| ^ ' 2,{l + w) 



IQ + Qw 

3(1 + u>) 



$(2) ^ ^(2) ^4^(1)2^^-1 






10H-6w 



3(1 -hw) 



a^^(i)a''^(i) 



3A^2^^^, 



a^(i)a*$(i) 



H±^5,vl;(i)a^vl;(i) 

3(1 + w) ^ 



(39) 

(40) 



In the case where the initial perturbations have been generated during a phase of one field inflation in slow-roll, it 
has been shown [Il| that, for the variable defined in Eq. ([55]) . T?.'-^-' ~ —27?.'^"'^^^ -I- O(slow-roll parameters). Under this 
hypothesis, we can use the constancy of 7?.'^' on super-Hubble scales and Eqs. pOMOp to derive the initial conditions 
satisfied by the two Bardeen potentials at second order at an initial time ?7init deep in the radiation era. Up to small 



corrections of the order of the slow-roU parameters, they are given by 

L J »7=i)i„it 

L J l)=»?init 

I 



a^.^(i)a*^(i) 



V=niait 



(41) 
(42) 



We also need to determine the initial conditions for the 
energy density contrasts and the velocities of the different 
matter components. Again, we assume adiabatic initial 
conditions, which means that the total fluid behaves like 
a single fluid. While at linear order the pressure and 
density perturbations are simply related by the sound 
speed 



5P(1) = c2<5p(i) 



at second order we have 



5P(2) ^ c2<5p(2) + i^^Sp) 



(43) 



(44) 



This implies that the adiabaticity conditions at second 
order reads 



where use of the first order adiabaticity condition was 
made to get the last equality. It can also be shown from 
the perturbation equations that the condition (|45|l re- 
mains valid on super-Hubble scales, hence giving a con- 
servation of the baryon-photon entropy on large scales, 
exactly as at linear order. 

Following the same procedure as at first order, we first 
use the Poisson equation (P^ and Eq. (P5|) at second 
order to get that 

SPimmt) = -2$(2)(^.^.^) + 8a>(i)(ryinit)$(''(^init) , 



and 



fci;^2)(k,77init) = -fc»7init 


_$(2)(k) 


+ <i>W(ki)$«(k2)f2-3^-f^-F^') 
V k \ki k2j 


Then using Eq. (|45|) we conclude that 


W?H?7init) = v'^'' (riinit) ^ vP {rjinit) , 


Si'Hmnit) = <5(2'(^i„it) 


3 

~ 4 


"^(2) _$(1)$(1)" 



JVinit 



(46) 



This completely fixes the initial conditions for the set of 
second order perturbation equations. 



E. Integrating the evolution equations 

Working in Fourier space, we first integrate the first 
order equations, that is 

• Eqs. (|17m8p for CDM and baryons assuming a 
source term of the form ((20)) for baryons; 

• the first eight moments of the Boltzmann equation 
for the radiation including the contribution of the 
polarisation. These equations are detailed in Ap- 
pendix El 

• the Einstein equations (|22ll25p . 

The initial conditions are detailed in iJ IIID 11 Technically, 
we have recast all these equations in order to use y as time 
variable and we remind that the time of decoupling is of 
the order of y = 3. We also remind that the modes of 
interest, that is fc > kcq., becomes sub-Hubble at y < 1. 
This first integration thus allows us determine all the 
first order perturbations as a function of time and wave- 
number. An example of the results of the first order 
integration is presented on Fig[7| f Appendix ICj) . 

Now, at second order, we integrate the same system of 
equations but supplemented by the source terms which 
are determined by the solutions of the previous integra- 
tion. We thus solve 

• Eqs. (jlTHlSp for CDM and baryon with the source 
terms (pS)) and (P7)) respectively. We recall our 
main hypothesis which states that the coupling of 
baryons to radiation can be described by the inter- 
action term ([50]) : 

• radiation is described by the first four moments of 
the Boltzmann equation with the hypothesis (|30p 
for the collision term; 

• the Einstein equations (|22j|25p with the source 
terms dSDMl)- 

The initial conditions are detailed in § HID 21 so that we 
are finally able to compute the evolution of the pertur- 
bation variables with y for any wave-number. 

Fig.[T]shows the evolution of the two second order grav- 
itational potentials. In particular, it shows that the so- 
lution is driven, as expected, toward ^^^^ = VP^^). It 
is to be noted that the convergence takes place before 
equivalence, at y = 1, as stressed in the following. On 
the other hand Fig. [2] depicts the evolution of the veloc- 
ities and density contrasts and shows that the photons- 
baryons plasma can safely be described as a single fluid, 
almost until the decoupling (shaded area on the figure). 




FIG. 1: Evolution of the two second order gravitational po- 
tentials, (;i'^'(ki,k2) andifj {i^i,]i.2) for ki — k2 — lOfcoq (top 
panel) or k^ = k2 — 20fcoq (bottom panel) and ki • k2 = 0. 



III. ANALYTICAL INSIGHT 

Before we proceed to describing the outcome of our 
numerical integrations, and in order to gain some insight 
into the physics of this intricate system, we present some 
analytic descriptions of its solutions. 



A. Heuristic argument and hypothesis 

Let us first assume that /b <C 1 so that the universe 
is mainly dominated by non-interacting cold dark matter 
and radiation components. When, in the radiation era, 
the CDM component is completely negligible the gravi- 
tational potential is determined by the density contrast 
of radiation. The latter however develops oscillations af- 
ter Hubble-radius crossing while those in the CDM fluid 
increases. It follows that, while still formally in the radia- 
tion era (pr > Pc), the cold dark matter component is ac- 
tually driving the gravitational potential. It then acts as 
an external driving term in the evolution equation of ra- 
diation. In such a scenario we then expect non-linearities 
that develop in the CDM sector to be transferred first to 
the gravitational potential and then to the radiation den- 
sity fluctations. 

To make this heuristic argument more quantitative we 
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FIG. 2: Top: Comparison of the baryons and photons velocity 
perturbation at order 2 for ki = k2 = lOfccq and ki ■ k2 = 0. 
It shows that Vr ~ i^^ with a good approximation until 
decoupling. Bottom: l.h.s and r.h.s of Eq. (|45p for the adia- 
baticity condition at order 2. It can be seen that this adia- 
baticity condition holds until recombination, hence justifying 
the approximation of ij lllll 



will thus assume that 

• we can first study the CDM-plasma system to de- 
termine the gravitational potential at second order 
where for simplicity the plasma is assumed to be 
radiation dominated; 

• then study the acoustic oscillations of the baryon- 
photon plasma driven by the gravitational potential 
derived this way, both at first and second order in 
the perturbations. 

For the sake of simplicity we will work in the tight cou- 
pling approximation. We recall that the time of decou- 
pling is both the time at which this tight coupling regime 
ceases to be valid and the time at which the radiation 
temperature is observed. 

The tight coupling approximation amounts to saying 
that the coupling terms !Fr and !Fh are so large that radi- 
ation and baryons behave as single fluid. It ensures that 
the two fluids have the same peculiar velocity {v^ = Wb) 
and implies that the anisotropic pressure of radiation 
vanishes {n^ — 0). 
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From Eq. p7|) . it implies that 



5b - l^r 



(47) 



At linear order, eliminating Ta in Eq. (|18p for radiation 
and baryons leads for the photons-baryons plasma to the 
continuity equation 

5;i + 3H(c2pi-u;pi)<5pi + (l + u;pi)(Ai;pi-3*') = , (48) 

and the Euler equation 

(? 
«;, + n{\ - 3c2^p,)«pi + $ + T^^^'^pi = , (49) 

where we have introduced the density contrast of the 
plasma b^\ = ^p^xj p^\ with 

(5/9pl = ipr + 5/9b , Ppl = Pr + Pb ■ (50) 

In the particular case at hand, it reduces to 

1 + i? 



5pi = 



l + |i? 



5r 



and the velocity perturbation are given by 

Wpl = t^r = Wb • 

The equation of state and sound speed of the plasma are 
easily obtained from the fact that P^\ = P^. They are 
explicitly given by 



1 



1 



Wpi 



3 + 4i? ' 



^s,pl 



3(1 + i?) 



(51) 



and are time-dependent quantities (simply because the 
relative contribution of the two components changes with 
time). 

At second order, the density contrast and velocity per- 
turbation of the plasma are given by 

(2) (.i+itjOr --Or (2) (2) (2) 



-"pi 



i + fi? 



X'pl' ^v-'= Vi, 



It can be checked that the plasma follows the equa- 
tion (|15|) and (|15)l but supplemented with the source 
terms 



'-'c,pl — *^c,a— pi '-'ciPl — *-'e,a— 



=pl, 



(52) 



where S'c,a=pi and S'e,a=pi stand for Sc,a and Se,a in which 
we take the fluid to be the plasma, that is (5a = (5pi, Wa = 
Wpi, etc. 

The validity of the tight coupling approximation at 
first and second order can be checked from our numeri- 
cal integration (which indeed does not make this assump- 
tion) respectively on Fig. [5] and Fig. [5] for the second and 
first orders. 



B. CDM-radiation system 

1. First order 

Deep in the radiation era, the gravitational potential is 
mainly determined by the radiation density contrast and 
decays on sub-Hubble scales. The contribution of matter 
is negligible in the Poisson equation and it follows that 
the potential is given by 



$(fc,77) = 3$(fc,77i„it) 






(53) 



where ji is a spherical Bessel function of order 1 and 
X = krj. The density contrast of radiation is given by 
S^ — —2$ on super-Hubble scales (see iJ IIID II) . 

Let us now turn to the evolution of the CDM fluid 
during the radiation era. In terms of the variable y the 
continuity and Euler equations lead to 



2 + 3y 
22/(1 + 2/) 



5c = Sq,{y) , 



(54) 



with a driving force determined by the gravitational po- 
tential 



5*(y) = 3$ + 



6 + 9y 



22/(1 + 2/) 



$ 



1 + y \k 



i-cq 



$ , (55) 



where a dot stands for a derivative with respect to y. 
The general solution of Eq. ([51)) is of the form 6c{k, if) = 
A-|-i?lnx + (5part where (5part is a particular solution given 

by 



y 5*(fc,77')77'ln(^^^d77', 



Opart — 



where 5$ is given by Eq. ([SS]) as long as 5pc <C 5pr- 
For y <C 1, the contribution of the particular solution 
is neghgible so that A ~ — f *&(?7init) and B ~ Q. The 
solution (I53p shows that <i> vary mainly when 77 ~ k~^ 
so that 5part ^ A + Bh\x. A and B can be obtained 
semi-analytically and are well approximated by A ~ 6 
and i? ~ -9 so that 5c ^ <i>(A:, 77init)(-4.5 -I- 91na::). 

This solution is valid as long as 5pc ^ 5pr in the Pois- 
son equation. However, on sub-Hubble scales 5-^ remains 
constant while, as we just saw, 5c grows logarithmically. 
Their contribution in the Poisson equation then become 
to be of the same order when y ~ y*{k), where y* is so- 
lution of 2/4-4-5 + 9ln{V2y^k/keq)] ~ 6/(1 - /J, where 
use has been made of a; = ^/2ky/kcq as long as y <C 1,. 
The solution of this equation is depicted on Fig. [3l For 
most of the scales of interest, i.e. for k ^ k^^., the contri- 
bution of the CDM in the Poisson equation is dominant 
before equality, i.e. y^ < 1). 

For these modes, which became sub-Hubble during the 
radiation era, we shall consider that CDM dominates in 
the Poisson equation and neglect the contribution of the 
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FIG. 3: The time at which the contribution of cold dark mat- 
ter and of radiation are comparable in the Poisson equation 
as a function of k/kcq. fceq is the wave- number of the mode 
that becomes sub-Hubble at the time of quality. For most of 
the scales of interest, i.e. for k ^ fccq., the contribution of 
CDM in the Poisson equation is dominant before equality, i.e. 

y* < i). 



radiation density perturbation, so that the Poisson equa- 
tion takes the form 



k^<,^-'^'-^^^ 






(56) 



Neglecting the contribution of baryons, since their den- 
sity contrast cannot grow because they are tightly cou- 
pled to the radiation, Eg. (jM)) then takes the form of the 
Meszaros equation [3^ 



Oc + :^^^ rfJc 



-5c = ■ (57) 



2yil + y) 2yiy + l) 

Its two solutions are a growing mode 

D+{y)^y + 2/3 (58) 

and the decaying mode 



D.{y)^-2^/{lTV) + D+iy)ln{^^±\] . (59) 



2. Second order 

At second order, as at first order, while deep in the ra- 
diation era the gravitational potential generated by the 
density contrast of radiation decays when a mode be- 
comes sub-Hubble. At this order the gravitational po- 
tential satisfies 



with 






Sr = S2- ^Si + ^ASs + HS'^ , 



(60) 



where the source terms are given by Eqs. pHI34p . Up to 
a fast decaying solution, the general solution is 



*(77)-3*(fc,77i„it) 






^init 



(61) 



with the Green function 
^1 



Gr{k,ri,ri') = — -I (cl^^xx' + l) sin[csj{x ~ x')] 

— Cs,rix — x')cOs[Cs.rix — x')] > . (62) 

On sub-Hubble scales, the leading terms in Sr are those 
quadratic in the first order velocity, of the form oc 
Ti^d'^vdiV ^ rj^"^, which behave as 77"^. All other terms 
in Sr behave, at best, as k~^rj~'^. Using the first order 
solution, the second order gravitational potential asymp- 
totically behaves as 



^ 



f{k,^) c. -:^c 



1 



81 *(l)(fcl,7yinit)*(lHfc2,?7init) 






ki • k2 



^2/ \ ^1^2 

X [cos {cs^-ckiT]) cos {cs^rk2'n) ~ COS (cs^i-kr]) 

-^sin(cs,rfci77)sin(cs,rfc2?7)] ; (63) 

and it can be checked that this term is indeed regular in 
/i^ = 1. Taking the homogeneous solution into account, 
^(^' decays as {kri)~^ on sub-Hubble scales. 

Now, the evolution of the density contrast of CDM fol- 
lows, using y as the time variable, the evolution equation 



J(2) 



3y 



■e 



-3^(2) 



22/(1 + 2/) 

6 + 9y 



22/(1 + y) 



^(2) 



1 



^cq 



(64) 



$^'^ = Ss^ , 



where the source term is given by 

Ss.^ 



1 



i-cq 



1 



2/ 



Sc 



Sc 

y 



1 



k 

y \kcq 



s. 



As in the previous section for the first order perturba- 
tions, CDM density perturbations grow faster than those 
of radiation so that, for any mode k that became sub- 
Hubble before ij^q, there exists a time of order yi,[k] such 
that for y > yi,[k] the gravitational potential at second 
order is determined by the cold dark matter. Where- 
after, even though we are still in the radiation era, we 
can neglect the contribution of the density perturbation 
of radiation so that the Poisson equation becomes 



fc2$(2) ^ ^2^(2) ^ _^(L^fc2 .(2) _ (g5^ 

42/ 



In this regime, the evolution of the density contrast of 
CDM at second order can be derived from a second order 
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Meszaros-like equation, in a similar way as at first order. 
Using Eq. ([55)1 and the fact that the main contributions 
in Ss^ in this regime come from 



Sc ^ -2d, {Sd'v) , d,Se ^ ~2 {djvd^d.v 
Eq. ((64|) takes the form 

A(2) , ^ + ^^ ki-i) _ 3(l--/b) rf2) 

^ ^22/(l+y) ^ 2y(l + y) ^ 

with 



5* 



M 



(66) 



(67) 



5.r = C< 



25c6c + 25l 



SrSr 



+ 



dcSc + 2Si + 



SrSr 



2y{i + y) 



yii + y) 



ki •k2 



,klk; 



2r2 



"'1"'2 



jL2 ' jL2 

'^l '^2 . 



Sim 



Now we shall neglect the effect of baryons, that is /b. 
This equation could then be called the second order 
Meszaros equation and describes a growth of CDM den- 
sity perturbation in a regime where radiation dominates 
the dynamics of the background while its density pertur- 
bations are negligible in the Poisson equation. 

The Green function associated to this equation is ob- 
tained to be 



3./ 



G(2/,y') = -2/'v/r+7(2 + 3y)(2 + 32/')x 



(69) 



yi + u 1 yi + M + i 

2 + 3u ^ 6 ^^ VTTH-l 
so that the general solution of Eq. (pT]) is 

6i'^ = CiD+iy) + C^D.iy) + f" G{y,y')SMiy')dy' . 

Jo 

In the limit where y ^ 1, y' ^ 1, the Green function 
behaves as 



Giy,y')^-y 



and the source term as 



l-(^ 



5/2' 



S{h)S{k2) 





jjlj — iv^ji yi^i, 


-. y2 ' 


with the Kernel 




i^(ki,k2) = 


'5 1, , 
7 + 2^^^ • ^^ 


fl 1\ 2(ki-k2)2" 



(70) 

In the limit y 3> 1, the particular solution dominates and 
our solution converges toward 

ij(2)(fc) ^Cif(ki,k2)5c(fcl)5c(fc2) , (71) 

that is toward the standard result (O describing the col- 
lapse of cold dark matter in a matter dominated era. 



f^L(ki,k2)0 




Logio[y] 



FIG. 4: (solid line) The second order potential computed in 
the tight coupling limit as a function of time y and of 6, angle 
between the wave vectors. It is compared to its expected late 
time behaviour (|72p . It is to be noted that the convergence 
toward this solution is extremely rapid and takes place as soon 
as equality is reached, e.g. y = 1. The results correspond to 



fci 



12A:cq. The difference in the amplitude 



of the function is due to the fact that the baryons component 
has been neglected in the derivation of Eq. I|72p 



The second order gravitational potential is then obtained 
from the Poisson equation 



i<I>(2)(fc,ry) 



-C^K(k,,k2) 



kik2T] 
k 



$(fcl)$(fc2), 



(72) 

up to terms of order 0{fh)- The convergence towards 
the solution ([72| is explicitly depicted on Fig. [4] where 
the behaviour of the exact (numerically integrated) sec- 
ond order potential as a function of time (and angle) is 
compared to its expected late time behaviour ([72]) . As 
detailed above, this solution is a better approximation 
for larger wave-numbers and at large y since we con- 
verge to this solution for y > y*{k). On this figure it 
can be observed that the convergence is extremely rapid 
and that the full kernel structure, including its angular 
dependence in ([7^ . is indeed to be observed in ^^^\ 



C. Baryons-radiation system 

We now want to understand the behaviour of the 
baryons-photons plasma, and in particular of its acous- 
tic oscillation, in the regime in which the gravitational 
potential is determined by the solutions of the previous 
section. 

We restrict our analysis to the tight coupling regime. 
And since it occurs for y < j/lsSj our solution will gain 
in accuracy when the period between CDM domination. 
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y = y^(fc), and the last scattering surface, y = yhss, is 
large, that is on the smallest scales. 



1. First order 

The computation at first order is well known '39i| and 
we review its main steps to compare to the more unex- 
plored second order case. 

In the fluid limit, the baryons and photons both obey 
a continuity and conservation equations p 711181 with a 
source term and in which the anisotropic stress of ra- 
diation can be neglected because of the tight-coupling 
approximation. As discussed in § IIII Al in the tight cou- 
pled regime, that is for k/r' ^ 1, it leads to the wave 
equation 



* 



R' 

'ITr 



$ 1 



* 



1 



* 



l + R 



(73) 



where R is defined in Eq. p3p , the sound speed is given by 
Eq. (|5ip . This is a wave equation with a forcing term on 
the r.h.s. which describes the oscillations of the plasma. 
For small wavelength modes, the variation of R and $ 
is small compared to the period of the wave so that we can 
construct an adiabatic solution by resorting on a WKB 
approximation; see e.g. Ref. ^391] for details. Defining 



©sw 
and the sound horizon 



^\5. + ^ 



1 



v/3[i + R{v')] 

the WKB solution takes the form 



dr/ 



Qswik^vi) 



(1 



cos[krs(ri)\- R<^ . (74) 



i?)l/4 

The velocity field, Ur = Wb can then be determined from 
the Euler equation ((49|) . 

This solution neglects the Silk damping effect that can 
be described by adding terms in fc^/r' in Eq. ([75]) so that 
the solution is exponentially suppressed by a factor 

I?(fc,7y) = exp ' 






where 



k-D^iv) - -, 



1 



l + R{i) 



16 
15 



R^jri') 

1 + R{-n') 



7V) 



so that 



e 



sw ■ 



i?$ ex exp 



/r2 



(75) 



15fc, 



eq- 



where the damping scale is of the order of kjj 

Since $ decreases as (fc/fcoq)~^: we conclude that for large 

wavenumber, Qsw{k,rii) goes rapidly to zero. 



2. Second order 

The former approach can be generalised at second or- 
der, but the behaviour of Qsw will change mainly be- 
cause the second order version of Eq ([75)1 has a r.h.s. 
which is steadily growing in the range of interest, i.e. 
after r^gq and on large k. 

At second order, Eq. (|73p will also contain terms com- 
ing from the second order Liouville equation [23| of the 
form 



5, 



pi 



ij ^ 



R' 

l + R' 



■Sr. 



fc2 



'e,pl ■ 



This source term involves terms which are quadratic in 
the fluid perturbation variables (i5, v) and the potentials 
($, ^). The former are exponentially suppressed due to 
Silk damping and the latter decrease as (fc/fccq)~^. We 
can thus neglect this source term as long as we focus on 
small scales. Defining, 



e(2) 



:>(2) 



SI- 



(2) 



.$(2) 



(76) 



the solution for Qgy^^ will be similar to Eq. ([75]) . When 
Silk damping is taken into account, 0^^^-|-i?$'^2) jg expo- 
nentially suppressed, exactly as at first order. The main 
difference with first order arises from the fact that the 
second order gravitational potentials are driven toward 
\Ii(2) ~ (j)(2) ^ (/c7y)2 after equivalence {y > j/eq) so that 
we expect that 



r(2) 



q(2) ^^^ , ^(2) 



^SW 



'R¥^^ 



(77) 



Now, the velocity of radiation is given by 



(1 + i?)w(2)J = -^ - (1 + i?)$(2) _ (78) 

Since 9^^-|-i?$(2) is exponentially suppressed due to the 
Silk damping, we expect that the Doppler contribution 
is also negligible. 

The equations (|72p and ([77)1 are the central results 
of the analytic insight of the non-linear regimes of the 
photons-baryons-CDM system at second order. They 
give the behaviour of the second order gravitational po- 
tentials at the time of decoupling together with the re- 
sponse of the photon-electron plasma. We also conclude 
that we expect Q^g^^ = -R^^"^^ to dominate the CMB 
temperature anisotropics on small angular scales (see Ap- 
pendix |B] for a discussion of the integrated Sachs- Wolfe 
contribution) . 

The bottom line of our analytic estimates is that on 
small angular scales, i.e. fc/fcoq ^ 1, the density per- 
turbation of CDM starts to dominate the Poisson equa- 
tion from y*{k) so that from this time to decoupling we 
can assume that the system is split in (1) the evolution 
of CDM and (2) the evolution of the photons-baryons 
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plasma which develops acoustic oscillations in the grav- 
itational potential determined by the CDM component. 
Because of Silk damping, Qsw + R^ dies out on small 
scales. At first order this implies that 6514/ ^ — i?$*^^^ 
which is suppressed by a factor k^^/k^ due to its evo- 
lution in the radiation era prior to yi,{k). At second 
order however, we still have for the same reason that 



e 



(2) 
SW„ 



_/j(j)(2) but now this term roughly grows as 
[krff' . Note that since fc^j is of order ISfcoq and that 
fceq roughly corresponds to a multipole I ^ 160, we ex- 
pect our analysis to give a good description of the system 
for £ > 2400. 



D. Comparison to numerics 

We now turn to the description of the numerical solu- 
tions of the system described in § III El Its solutions will 
be described in the light of the analytic description we 
just developed. 

Fig. [5] shows the result of numerical integrations for 
second order quantities of interest. They are compared 
to our approximate formula that, we recall, is expected 
to be valid in the tight coupling regime. It shows indeed 
that modes with k > fceq relax temporarily toward the 
solution ([77)1 . The exact solution exhibits though large 
oscillations that are thought to be due to the acoustic 
oscillations that are present in the plasma at first order, 
but their average turns out to coincide with the proposed 
analytic formula (long dashed lines) as long as a strong 
coupling is ensured. The Silk damping effect is observed 
to play a key role to actually damping the oscillations. 
This effect is all the more important that k is large. Note 
that the impact of the oscillations on the observational 
quantities will also be damped by the finite width of the 
last scattering surface. The wave-number corresponding 
to this width is of order of lOfccq which is smaller than 
the damping scale. 

When the coupling becomes loose, log]^o(y) approach- 
ing 0.6, the numerical solution departs from the expected 
solution (and it converges toward 0) as the full Boltz- 
mann hierarchy is now at play. We also depict the ve- 
locity term which can be checked to be negligible, as 
expected. 

From these set of results we can then argue that 
the approximate analytic solution described in (j77p cap- 
tures the physics of the dominant terms of the CMB 
anisotropics on small angular scales. Here we have explic- 
itly checked that this form is consistent with the physics 
of recombination when the collision effects are taken into 
account. We limit though the collision effects to their 
first order expression. We expect nonetheless that an 
exact calculation, up to second order, would not signifi- 
cantly alter our conclusions, the collision physics playing 
a role only during a limited period of time. Although 
this is certainly desirable to do such a calculation, this is 
beyond the scope of this paper whose goal is to estimate 
the order of magnitude of the non-Gaussianity on these 




FIG. 5: Behavior of ef^r (black) and kv^'^ /VS (red). The 
numerical integration is depicted with solid lines while the 
analytical estimates is plotted in dashed lines. From top to 
bottom, we have ki — k2 = 30 fceq. and ki — k2 = 40 fccq. . The 
vertical grey zone represents the "surface" of last scattering. 



scales. 

In the following we explore the observational conse- 
quence of such a finding on the temperature bispectrum 
at small scale. 



IV. SIGNATURE IN THE COSMIC 
MICROWAVE BACKGROUND 

A. Flat sky approximation 

Since our approximations hold on small angular scales, 
it is amply sufficient to treat the sky as flat to compute 
the properties of the CMB anisotropics. We thus decom- 
pose the CMB temperature anisotropics in 2D-Fourier 
space as 



e(n)= / |^e(i)e^ 



iln 



so that 



e(i) = /^e(n)e- 



iln 



(79) 



(80) 
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On the other hand 8(n) can be expanded in Fourier 
modes as 



e(n) = 



dk 



(27r)3/2 



0(k,77LSs)e- 



i(fer')LSS+£'LSski-n) 



(81) 



where I?lss is the angular distance of the last scattering 
surface given by -Dlss = J^" l/{H[u]u)du. \i.± is the 
projection of k on the sky, i.e. k = fc,.e + k^ where e is 
the direction of the (flat) sky. 

In Eq. ((8T|) we refer to 9(k, 77lss) as 



0(k,?/LSs) = / d77e(k,?7)i;(?7,ryLss) 



(82) 



considering the observed CMB anisotropies as a super- 
position of spheres of temperature anisotropy weighted 
by the visibility function f(?7,7yLss) which peaks at rj^ss- 
The angular power spectrum is nothing but the two di- 
mensional power spectrum of 0(1), 



(9(1)9(1'))= .5(2) (1 + 1')Q 



(83) 



We now need to determine 0(k, 77) in terms of the per- 
turbation variables. The CMB temperature anisotropies 
are usually split as an intrinsic Sachs- Wolfe effect, a 
Doppler effect and an integrated Sachs- Wolfe contribu- 
tion. As discussed in Appendix |B] the integrated Sachs- 
Wolfe contribution is expected to be negligible on the 
scales of interest in our study. 

At flrst order, the Fourier component of temperature 
anisotropy for a mode k, in a direction e emitted at a 
comoving distance 770— ?? is dominated by the Sachs- Wolfe 
and Doppler terms. 



9W(k,,y) = 9^^^(k,,y)-f9[,^^p(k,,y) 



(84) 



The Sachs- Wolfe term is related to the perturbation vari- 
ables by 

9lj^),(k,,y) = ^(k,,7)-fci>W(k,,7), (85) 

while the Doppler term is given 9dop = — wk.e so that 

e'i^liKrj) - iKvi'Hk,rj) , (86) 

with fcr = k • e. Using Eq. ((8T|) . we deduce thus that 



9(i)(l) 



2^^£ss 



dfc^5(i)(k)e''='-''Lss , (87) 



with 



5(1) (k)= /d77t;(77)g(i)(k,77)$(i)(k,0) 



(88) 



and where k_L = 1/Dlss- Using the definition of the ini- 
tial power spectrum ($(k,0)$(k',0)) = S'^^^k + k')P{k), 
we finally get that Ci is given by 



C, 



1 



2^^Ess 



dkrP 




g^'\k) . (89) 



This reproduces the main features of the CMB angular 
power spectrum, as checked on Fig. [HI 

At second order, and for the scales of interest (fc ^ keq 
and k > ko), we stress again that the knowledge of the 
exact expression of the terms quadratic in the first or- 
der variables in the second order Sachs- Wolfe effect are 
not needed since they are suppressed because of the Silk 
damping or because of the decaying of the potential dur- 
ing the radiation era. The analysis of § IIIIC 21 shows 
that the Doppler term is much smaller than the intrinsic 
Sachs Wolfe term of Eq. ([TT]) . Thus, the main contribu- 
tion to the second order temperature anisotropy is well 
approximated by 

9(2)(k,,y) ^ Qf^{k,v). 

C, ^ll^ + <i>(2)(k,,y)^_i?$(2).(90) 



We deduce that 



9(2) (1) = \^ /dM^'Hk)e"='-''-^ , (91) 



with 



5(2) (k)= f dr]v{r])e^^\k,T]). 



(92) 



It can be rewritten in terms of the initial first order grav- 
itational potential as 

g(^'Hk) = 2C 4f)(ki,k2)<l>W(ki,0)<I>(i)(k2,0) , (93) 
hence defining Z^^-*. 

B. Bispectrum 

In the fiat sky approximation, the reduced bispectrum 
bi^i^i^ is defined from the 3-point function as 

(9(li)9(l2)9(l3)) = {27r)-^6(^^\h+h + h)biM, ; (94) 

see e.g. Ref. [iJI- With the previous definitions, it can 
be expressed as 
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'■>hhh 



2t^^ss 



dkridkj.2 



/(f)(-ki,-k2)g«(ki)5«(k2)P 




(li ^ I2 ^ I3 ^ li) + (li ^13^12^ li) . 

I 



(95) 



C. Numerical computation 

In order to perform the previous integrals, we need to 
specify the initial power spectrum. We assume that the 
power spectrum is scale invariant and we normalize it 
using the results of WMAP, that is 



with 



P,„ . 2,''Pi (%) i- 



Ai = 3.33 X 10" 



(96) 



(97) 



We then compute the bispectrum of an equilateral config- 
uration for which all momcntums are equals; h = I2 — Is- 
The only free parameter for such configuration is the 
norm I of the three vectors. 
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FIG. 6: Black: the bispectrum for the equilateral configura- 
tion computed in the flat sky limit. The thin colored lines 
represents the bispectrum that would be obtained by assum- 
ing a constant initial f^^ and a linear transfer function, that 
is neglecting the non-linear dynamics. From bottom to top 
we have plotted /*^ = 1, 10, 100, 1000, 10000. 



The order of magnitude of the amplitude of the bispec- 
trum can be understood from the following rule of thumb 
for modes larger than kz). According to our analysis, con- 
sidering the equilateral configuration where fci = ^2 = fc, 
the second order temperature anisotropy on the last scat- 
tering surface is of order 



1 



e(2)(k,?7LSs) ^ -^i?LSS$<'^(k,r;LSs) 



-Rlss 



2yLss_fc^ 



<i>'-^\k,0)TJ^\k) 



where we have assumed that, on average, the Kernel 
ii'(ki,k2) is of order unity. Now, assuming a constant 
primordial /^j^ evolved with the linear transfer function, 
the second order temperature anisotropy would roughly 
be of order 



ie(2)(k,,7Lss) 



i?LSs/NL[<i>('nfc,0)l T«(fc) (98) 



since, for these modes, the integral on the visibility func- 
tion keeps only the average of the Sachs- Wolfe contri- 
bution. The ratio of the contribution of the non-linear 
dynamics compared to a primordial non-Gaussianity is 



2yLSS 



^eq 



r 
■'0 



(1) 



(fc). 



(99) 



Since for large modes the gravitational potential has been 
decaying as {kT})~'^ in the radiation dominated era, and 
growing logarithmically when the potential started to be 
determined by the cold dark matter component (see the 
analysis of section UlIB ip . we deduce that the first order 
transfer function is typically given by 



T^'\k)^A{k)[^-f 



(100) 



where A{k) is a steadily growing function. At the Silk 
damping scale we find numerically A{k£i) ~ 10. We thus 
conclude that in the bispectrum, the evolution for I ~ 
koDi^ss — 2400 is equivalent to a primordial 



The result is depicted on Fig. [5] and is compared to the 
bispectrum one would obtain from a initial constant f^^ 
assuming a linear transfer function. It appears that on 
scales that range from I = 1000 to I = 3000 the bispec- 
trum resembles that of an effective constant primordial 
/* of order 25. 



/»L - ^VLSsAiko) ^ 25 

evolved linearly. This estimates of the order of magni- 
tude is in complete agreement with Fig. [S] where it can 
be read that for multipole ranging from 2000 to 3000 
the amplitude of the bispectrum is comparable with the 
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one that would be obtained from a constant f^^^ ranging 
between 10 and 50. 

There is no guarantee however that for arbitrary ge- 
ometries the shape dependence of the temperature bis- 
pectrum would be that of a constant /^^j^. It is rather 
determined by the kernel shape of the form ([77|) . 



V. CONCLUSION 

This article investigates the non-Gaussianity that 
arises in the CMB temperature anisotropics due to 
the post-inflationary non-linear dynamics during the 
radiation and matter dominated era. More specifically it 
aims at identifying the leading mechanisms and leading 
terms that determine the shape of the CMB bispectrum 
on small angular scales. 

The driving idea that we have pursued throughout the 
paper is that at small angular scales the second order 
CMB anisotropics trace the second order gravitational 
potential as it is shaped by the CDM component during 
its sub-Hubble evolution. To give support to this picture, 
we have developed both analytical insights into the joint 
evolution of the density potentials and the temperature 
fluctuations and numerical tools. 

We have thus solved numerically the joint evolution 
equations of the cosmic fluids up to second order. We 
have been able to check that at the time of decoupling, 
the second order potential indeed traces its expected 
shape. This conclusion is summarised and illustrated 
on Fig. m At this stage, and as long as one restricts 
these results to the tight coupling regime, no approxi- 
mations have been made. The accuracy with which the 
fc-dependence of the matter dominated mode coupling 
kernel, i.e. Eq. ([6]), is recovered is truly remarkable. We 
stress that this is due to the fact that for the physics at 
work at small scales, i.e. k/k^q S> 1, the density per- 
turbations of the CDM component start to dominate the 
Poisson equation much before equality. This implies that 
the non-linearities developed by the CDM can be trans- 
ferred very efficiently to the gravitational potential even 
before the beginning of the matter era. 

Determining exactly how this mode coupling ker- 
nel is actually transferred to the source term of the 
CMB anisotropics relies on further numerical integrations 
through the recombination era. At this stage, the only 
approximation we make concerns the Compton scatter- 
ing collision term entering the Boltzmann equation for 
radiation at second order. We did not use its full second 
order expression but we argue that it can be reduced to 
its formal first order form, namely to Eq. (|30p . Such an 
assumption is clearly valid in the tight coupling regime. 
Actually this is the only term appearing in the collision 
term in the baryons rest-frame as long as tight coupling 
at first order is efficient. We then argue that when the 
coupling drops. Silk damping effects effectively suppress 
all other contributions. 



This leads to the behaviour depicted on Fig. [5] for the 
main source term of the temperature anisotropics. In 
particular 9^^, the monopole of the second order source 
term, is found to be attracted toward a non-vanishing 
and non-oscillatory term, —R^^"^', where we recall that 
R is the baryon to radiation ratio, while the dipole con- 
tribution, and thus the Doppler effect, vanishes. As a 
result the main contribution to the CMB temperature 
anisotropics at second order is found to be directly pro- 
portional to the second order gravitational potential. It 
has to be noted that the efficiency with which the sec- 
ond order term converges to this form is considerably 
accelerated by the Silk damping effects which efficiently 
suppress the oscillatory parts of the solution. This results 
are clearly illustrated on Fig. [51 We argued from order of 
magnitude arguments that, because the damping scale is 
typically of order 15fccq, our description shall be valid for 
I > 2400. This observation is the basis of the main result 
of this paper. Actually,as Fig. [5] tends to show, it seems 
that this description may be valid at lower multipoles. 

Finally we explore the consequence on the CMB 
bispectrum. For obvious reasons we use the small angle 
approximation to perform the numerical integrations. 
The bispectrum for equilateral configurations is illus- 
trated on Fig. [B] Wc show that for these configurations, 
its amplitude corresponds to what a primordial non- 
Gaussian potential of /* of order 25 would have given 
(also for an equilateral configuration). As shown in the 
text, this number can easily be recovered from back-of- 
the-envelop calculations. The first lesson that can be 
drawn from this result is that it gives a signal larger than 
what a model with a primordial /^^ of order unity would 
give! The second lesson is that the /-dependence of the 
bispectrum is expected to be different from the one 
induced by primordial mode couplings. It is expected 
to have a specific shape as encoded in the CDM kernel 
expression. 

In conclusion, this work offers a breakthrough insight 
into the physics of CMB in the non-linear regime and on 
small angular scales. It identifies what is, as we argued, 
the main small scale contribution of the bispectrum, 
hence filling the gap with the standard results that 
have been obtained in the weakly non-linear regime 
of gravitational clustering of dark matter. We did 
not check this result against a (yet non-existing) full 
second order Boltzmann code, and this is probably 
desirable, but we argue that, given the amplitude of 
the effects, all other contributions will be subdominant. 
With such a large signal, detection of this bispectrum 
should be easily within reach of future CMB experiments! 
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APPENDIX A: DESCRIPTION OF RADIATION 



To describe the evolution of radiation, we use the first 
moments of the Boltzmann hierarchy (see e.g. Refs. |44| ) 
including polarisation. The hierarchy reads 



e', = k 



E' 



i 



2i- 1 
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1 



2^ + 3 



e 



e+i 



e^-fe— (e2-\/6-B2 
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-Eg-i 



Vi^+W 



21-1 2i + 3 
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(Al) 



(A2) 



where the first moments are related to the fluid variables 
by 



eo = i<5. 



81 = —kv-c 



62 = —k^TT, 

12 



The Boltzmann hierarchy is infinite and wc truncated it 
after the multipole £ = 8, when computing the first order 
and after d — Z when computing the second order. In 
order to cut the hierarchy without numeric reflection j45|] , 
we use the free-streaming solution of this hierarchy, and 
use it to express in the last equation the multipole (. + 1 
in function of the multipoles for I and I — 1. Explicitly, 
the closure relation reads 



Ef,+i 



2^ + 3 



e, 



2£ + 3 



e, 



(A3) 



k-q 2i-l 

{2i + 3)VT+T ^ (2^ + 3) 7(7+ 3) (^ + 2) ^ 
fcTyVT^T ^ {2i-l)^{i-l){l-2) ^"' 



APPENDIX B: THE INTEGRATED 

SACHS- WOLFE EFFECT CONTRIBUTION TO 

THE SMALL SCALE BISPECTRUM 

At linear order, the contribution of the integrated 
Sachs- Wolfe effect on small scales is usually small because 
the time dependence of the potential vanishes in the mat- 
ter dominated era. This is no more the case at second 
order. It is thus legitimate to investigate the impact of 
the time dependence of the second order gravitational 
potential on the amplitude of the bispectrum. 

At linear order the expression of the temperature 
anisotropics is 



e(i)(i) 



V27r 

^LSS 



g(k) dkr 



(Bl) 



while at second order an extra source term should be 
included. It is formally given by 

©isw(l) = y d77-{^I'(2)[x(r7),r;]+<i>(2)[x(r;),,7]}(B2) 
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assuming instantaneous recombination at 77 = t^lss- To 
estimate the magnitude of this effect, we assume that 
v]/(2) — $(2) g^yj^ ^jjg^^ ^]^q tiffic dcpenclence is the one 
otained in Eq. ((7^ . 



where the spectral index Ug varies a priori from 1 (at very 
large scale) to —3 at very small ones - it is a priori of 
the order of say —2 at the scales of interest - then the 
integral over fc^ — fc^a leads to the factor, 



^LSS 



2 -'^-/v 1^ v^*-i5 ^*-z 1 '/ijDGy 1 X^'^J 



where the time dependence has been explicited. Obvi- 
ously, expression (jB2[) gives an extra term contributing 
to the bispectrum. Consistently with our former analy- 
sis, let us evaluate this contribution in the small angle 
approximation. This leads to 

/oo 
dfc,,dfc,,P(fci)P(fc2)XArL(ki,k2,77Lss) 
-00 



rvi.ss ^^ ^ 

X / d7jiv{'ni)dr]2v{T]2)^—g{ki,r]i)g{k2,7j2) 

Jo ^LSS 

X exp [i/cri?7LSS + ifcr2'7LSs] + sym. 



(B4) 



that has to be compared with Eq. (|95)) . The integral over 
77 can then be performed to give 



»)LSS 



X / dT]iv{r]i)dri2v{r]2)g{'ki,rii)g{'k2,ri2) 

Jo 

xw(fcri -I- fcr2,'7LSs) +sym. (B5) 

with 

/, ^ 1 - ifc?/LSS - exp(ifc77Lss) fr,^. 

w(k,rii^ss) = r^-^ • (B6) 

'^ ^LSS 

If one examines the UV convergence properties of this 
expression (for the integrals over fc^), it appears that the 
integral over fc^ -I- kr^ converges at a scale given by the 
inverse of t^lss, e.g. 

/oo 
dA:w(fc,77Lss) = ^ , (B7) 

-00 jyLss 

due to the oscillatory behaviour of w, whereas the inte- 
gral over fc^i — kr2 converges because of the power spec- 
trum shape and therefore at a scale which is of the order 
of /i/?7lss or /2/'7lss (whichever is smaller). 

Thus if the power spectrum is approximated by a 
power law. 



dkPii/^+pJ^s) = 

;2(--3)-i0F,^--r(|-n.) 



(B9) 



r(4 — Hs) 

that is 637r/256 (jyLss/O" for ?i, = -2. 

This is to be compared with the amplitude of the in- 
trinsic effects we have computed. The latter differs in 
Eq. ([55]) because of the absence of filtering function w. 
The amplitude of the bispectrum is then roughly given 
by, 

/oo 
dkr.dkr, P.S>{kl) P<S,{k2) 
-00 

xKNL{ki,-k2,VLSs) (BIO) 

SO that its amplitude is dominated by the square of 



dkr.dkr, P{ki) P{k2) Knl{^iM, ^LSs) / 

-00 J — ( 



:s"=;2(^-i)-i0^r(|-^) , ^ 

^,„ „_^ ^' — ,(B11) 



'^LSS' 



V{2~^] 



F<E,(fc)-fc" 



(B8) 



which is equal to {iir/Sf (77lss/0^° for n^ ^ -2. The 
ratio of the two contributions scales then as 1/{RI) in 
favor of the intrinsic effect. 



APPENDIX C: CHECK OF THE NUMERICAL 
INTEGRATION 

We report in this appendix the results of the first order 
numerical integration. We first report in Figs. [7] and [5] 
the evolution of the perturbed quantities where it can 
be seen that for y > y^(fc) the gravitational potential 
tends to be determined by the cold dark matter density 
perturbation. We also report the angular power spec- 
trum obtained from the flat sky approximation using the 
expression ([5^ . The linear dynamics is then used to cal- 
culate the bispectrum arising from a constant primordial 
f^^ evolved linearly. It can be checked that the form 
obtained on Fig. [H] is completely consistent with the lit- 
erature (see e.g. Ref. p6|). 
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FIG. 7: From top to bottom, left to right: evolution of the Bardeen potential, $, and the density contrast 5 for respectively 
radiation, baryons and cold dark matter. The solid line corresponds to fc = lOfcoq and the dashed line corresponds to fc = 20fccq. 




-0.5 0.0 

Logjo[y] 




FIG. 8: Left: Comparison of the baryons and photons velocity perturbation at first order for k — lOfccq. It shows that 



„(i) 



«Jj with a good approximation until decoupling. Right: Comparison of j5i and -^Sr ■ It can be seen that the 



1a(i) 



lx(l) 



adiabaticity condition holds until recombination, hence justifying the approximation of § IIIII 
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FIG. 9: Left: in solid red line, the angular power spectrum using our code and the flat sky approximation (which does not 
includes the late ISW). The black dashed line represents the spectrum obtained using CAMB, which also takes into account 
the late ISW, but no reionisation. Middle: The precision on small scale depends highly on the computation of the visibility 
function. On small scales, our transfer function is approximately 30% smaller than the one predicted by CAMB and is thus a 
good approximation. Right: the bispectrum obtained from a primordial constant f-^^^ — 100 evolved linearly. 



